Discrete Coordinate Scattering Approach to Nano-optical Resonance 
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In this letter, we investigate the coherent tunneling process of photons between a defected circular 
resonator and a waveguide based on the recently developed discrete coordinate scattering methods 
(L. Zhou et a/., Phys. Rev. Lett. 101, 100501 (2008)). We show the detailed microscopic mechanism 
of the tunneling and present a simple model for defect coupling in the resonator. The Finite- 
Difference Time-Domain(FDTD) numerical results is explored to illustrate the analysis results. 

PACS numbers: 



Recently, resonant optical cavities has attracted lots 
of interesting from both experimental and theoretical as- 
pect [1]. Devices based on these micro-cavities has been 
widely used for processing optical signals 0. In most 
hardware of these signal processing, the elementary unit 
consists of a micro-cavity and a side coupling waveg- 
uide [l| . This unit has been experimentally realized based 
on photonic crystal, most recently ring made with GeN 
nanowire Q. With such kind photonic setups, photon 
blockade was experimentally observed in this unit [J]. 
For this unit, it is well known that the incident wave is 
totally reflected by the ring resonator when it is resonant 
with the modes of the ring. Thus, the resonator takes 
the role of frequency-selection [Ej] in this unit. 

Such transfer process in the unit is crucial to design 
more complicated devices. To investigate this phenom- 
ena, one way is to directly solve the Maxwell equa- 
tions. However, it is impossible to get an analytic re- 
sult to reveal the key physical mechanism for even a 
simple configuration. Thus, this is usually done by 
numerical simulation, such as Finite-Differential-Time- 
Domain(FDTD) method @, which is used as a validation 
approach. Theoretically, coupled mode theory(CMT) is 
introduced to study the resonator scattering properties 
in this unit [H, Q- However, this is done in the mode 
space, and therefore not intuitive in real space. 

In this letter, we generalize the most recently devel- 
oped discrete coordinate scattering methods @ to inves- 
tigate the detailed scattering properties for the optical 
signal processing. Here we model the resonator and the 
waveguide as array of coupled cavities. And the coupling 
between resonator and waveguide is modeled as hopping 
of the photon among cavities. We present the method by 
a qualitatively analysis of an example of a notched-ring 
resonator. Rigorous FDTD numerical simulations is also 
performed to verify the obtained analytical result. 

The elementary units under investigation is schemati- 
cally illustrated in Fig. Ufa). The defected ring resonator 
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FIG. 1: (Color online) Frequency-selective filter. (a) A 
notched ring is evanescent ly coupled to the waveguide. The 
notch couples the CW and CCW traveling modes; (b) Model 
of the unit by discrete coupled resonator array. The notch 
is characterized as point inhomogeneous of the tunneling 
strength along the ring. 



is evanescent ly coupled to the signal waveguide. For 
perfect ring, it supports the clockwise(CW) and counter 
clockwise (CCW) traveling wave and there is no interac- 
tion between the two traveling wave modes. However, 
the scattering of the defected notch causes the recom- 
bination of the CW and CCW modes. In the context 
of CMT [5[, the notched ring is mapped to the coupled 
double-ring resonator and the notch in single ring plays 
the role of the gap, which couples the double rings. Here, 
we treat the defected circular resonator as the discrete 
coupled cavity array with coupling constant different at 
the defected area, which is illustrated in Fig. Hfb). The 
array consists of N resonators with tight-binding cou- 
pling. The notch in the ring is marked with a dashed 
line. The Hamiltonian here reads 



J ^ajaj+i + h.c.^ + S J ^ajai + h.c. 



* Electronic address: dhui@itp.ac.cn 
t Electronic address: suncp@itp.ac.cn 



(i) 

where aj is the creation operator for photon with en- 
ergy uo the j — th cavity and J is the tunneling strength 
between nearby cavities with defected induced inhomo- 
geneous coupling SJ. By the Fourier transformation 
Oj = ^2 exp (— ijk) c\,/\ /r N, the Hamiltonian in the mo- 
mentum space reads H c = (u — 2 J cos k) c\ck + 
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FIG. 2: Ez-field amplitude of ak(B\ue line) and /3fc(Red 
dashed line) mode on the notched ring. The position of notch 
is marked with dot line. For mode, the amplitude of Ez- 
field at the notch is at maximum. However, the amplitude is 
null for /3k mode. 



5 J J2k p ex P (}P) c k c P + ^' c - /N- In the rotation wave 
approximation, we keep only p = ±fc term and obtain 
the Hamiltonian as H c = ^Z k>0 £k (c\ c k + ct fc c_/^ + 

^2 k>0 25 J (exp (—ik) c\c- k + h.c.^ /N, where s k = co - 

2 (J — SJ/N) cos/c. Physically, c^and c_ k stand for the 
CW and CCW traveling wave, which has been introduced 
phenomenally in CMT. The coupling between the two 
modes is induced by the defect of the inhomogeneous 
of tunneling 5J . Diagonalizing the Hamiltonian with 
unitary transformation cJ = exp (ik/2)(al + pl)V2 and 
c -k = ex P (—ik/2) (a\ — Pl)V2 for k > 0, we obtain 



H, 



/c>0 



where e k = (sk ± 25 J/N) and 

a i = 4=E 2fl i co8 Kj'- 1 /2)i] 



i=^2m]sin[(j-l/2)fc]. 



(2) 

(3) 
(4) 



are the eigen-modes of the notched ring. And they were 
related to photonic molecule in Ref. 0, [l2[. Actually, 
this single defected notch induced the split of the degen- 
erate resonant peak of CW and CCW traveling wave. 
The splitting can be probed by checking the transmis- 
sion rate in the waveguide, which will be investigated in 
the following parts. Actually, those modes correspond 
to a symmetric and antisymmetric state with respect to 
the notch. The amplitude for the modes otk and /3k is 
illustrated in Fig. [2j In this discrete model of the ring 
resonator, the position jnotch of notch should between the 
resonator ao and ai, namely jnotch = 1/2. Thus for the a 
mode, the amplitude of the field at the notch is at max- 
imum. However, there is no field in the notch for the f3k 
mode. 

To probe the optical property of the resonator, a trans- 
mission line waveguide is introduced to couple the signal 
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FIG. 3: Plot of the reflection rate with the parameters uo — 
U! = 8.772 x 10 14 Hz, J = J' = 3.3 x 10 14 Hz, g = 3x 10 13 Hz 
and 5 J — 2.98 x 10 13 Hz. (a) Theoretical calculation of the 
reflection of the filter; (b) FDTD simulation of the reflection 
and transition rate. Two resonator peaks appear at the cavity 
mode ei . 



in the microcavity. Physically, the waveguide can be opti- 
cal fiber, nanowire and defected lines on photonic crystal 
for different systems. Here, we also model the trans- 
mission line waveguide as the discrete coupled resonator, 
which is illustrated in Fig. Hfb). Each resonator couples 
to the near-by resonators by means of photon hopping. 
Those coupled-resonator array has been proposed to be 
a new type of optical waveguide [9j. The Hamiltonian for 
the transmission line waveguide here reads 



#z = -J'^(^ i+ i+h.c.) 



■u'b)b h 



(5) 



where is the bosonic creation operator for photon on 
the j — th resonator. The tight-binding Hamiltonian here 
is a good approximation for modeling the coupled res- 
onator array, which is proved in Ref. [9]. And in Ref. 
[8[ , it is showed the model is equivalent to the simple di- 
electric waveguide under the linear dispersion area with 
right going wave 4)r(x) and left going wave 4>l{x) as 
showed in Ref. [10[. The waveguide couples to the cav- 
ity at the 0-th resonator as the tunneling between 0-th 
resonator and the jo-th resonator of the ring resonator 

Hj = g (^Oj Q bo + aj bl^J. The total Hamiltonian for this 

elementary unit reads H = H c -\- Hi + Hj. For the single 
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photon case, the eigenwave stationary scattering function 
can be written as 

j k>0 

where Uj is the amplitude of photon in the j — th res- 



onator on the waveguide and Vj(wj) is the amplitude 
of photon on the resonator with mode c^(ft). Substi- 
tuting the stationary scattering function into the time- 
independent Schrodinger equation = E\<&), one 
can obtain equations for the amplitude of the site as 



J' (u j+1 + uj-!) + Juj + A> E K + «*) e-^/a-io)* + („ fc - «, fc ) e-**Oo-i/2) 

v 2j?v fe>0 



e + Ufc + f e <(V2-io)* + e iOo-l/2)A 

V2iV V / 

£ ~ Wk + ( e i(l/2-io)k _ e iUo-l/2)k\ 



Euj, 

Ev k , 
Ew k . 



(7) 

(8) 
(9) 



When the incident energy coincides with the energy of 
a mode, namely E = e^, the amplitude i/o = by Eq. dHJ 
and k;^ = by Eq. ([9]). Thus only a mode is excited here 
and the Ez-field at the notch is at maximum. When the 
incident energy is resonant with the energy of (3 mode, 
namely E = e~j~ , the amplitude uq = by Equation® 
and Vk = by Equation®. Here, only (3 mode is excited 
and the amplitude of the field is zero. This phenomena 
has been numerically predicted in Ref. [Bj]. 

For the scattering process, the photon incidents 
from the left side of the line waveguide and is par- 
tially reflected by the notched ring resonator ar- 
ray. The amplitude in the waveguide reads Uj = 
[exp (iqj) + r exp (-iqj)} 6 (-j) + s exp (iqj) (j) , where 
r(s) is the reflection (transmission) coefficient of the pho- 
ton. The coefficient of the transition here is obtained 

K 1 "^)" 1 ' (10) 

where Q = 2g 2 /N £ fc>0 cos 2 (£(j - 1/2))/ (E - e+) + 
sin 2 O(j - 1/2))/ (E - e~). The reflection rate 1 - \s\ 2 
peaks at the resonance with the internal ring mode sj^, 
where the optical signal is totally reflected. To see this 
resonant filtering effect, we plot the transmission rate |s| 2 
verse the the incident energy of photon in Fig. [3](a). For 
comparison, a rigorous FDTD[6[ simulation is performed 
with free available code F2P pj| for Transverse-Magnetic 



mode. The simulation is performed with the following 
parameter: no = 1, n w = 3.2, R\ = 5a, R2 = 4.4a, 
= 54° and SO = 1.5°, where a = 0.25/ira. The 
transmission and reflection rate is plotted in Fig. 02(b) 
around the resonant frequency Sk = 0.331 x 2irc/a. 
In Fig. [3](a), the parameters for this defect ring res- 
onator unit are chosen to fit the numerical simulation: 
co = u) X = 8.772 x 10 14 i^, J = J' = 3.3 x 10 14 i^, 
g = 3 x 10 13 Hz and SJ = 2.98 x 10 13 Hz. Therefore, our 
model confirms the function of the filter and its mech- 
anism. However, we should notice a drawback to this 
approach: the quantitative analysis will relies on the nu- 
merical simulation. At present, the coupling constant 
at the notch and between the transmission line and the 
resonator should be fitted by numeric simulation. 

In conclusion, we have presented a theoretical inves- 
tigation of the notched ring filter with discrete coordi- 
nate scattering methods. The detailed mechanism of the 
coupling between the clockwise and count-clockwise is 
revealed to be the inhomogeneous tunneling in the res- 
onator array system. We have predicted some experi- 
mentally accessible results for the coherent transmission 
and reflection property and present detailed explanation 
to field amplitude distribution at the notch. 
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